Hysteresis and noise in ferromagnetic materials with parallel domain walls 
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We investigate dynamic hysteresis and Barkhausen noise in ferromagnetic materials with a huge 
number of parallel and rigid Bloch domain walls. Considering a disordered ferromagnetic system 
with strong in-plane uniaxial anisotropy and in-plane magnetization driven by an external magnetic 
field, we calculate the equations of motion for a set of coupled domain walls, considering the effects 
of the long-range dipolar interactions and disorder. We derive analytically an expression for the 
magnetic susceptivity, related to the effective demagnetizing factor, and show that it has a logarith- 
mic dependence on the number of domains. Next, we simulate the equations of motion and study 
the effect of the external field frequency and the disorder on the hysteresis and noise properties. 
The dynamic hysteresis is very well explained by means of the loss separation theory. 



PACS numbers: 75.60.Ch, 75.60.Ej 



I. INTRODUCTION 



The study of ferromagnetic hysteresis represents an 
open field of current interest, due to the applications in 
magnetic recording technology and spintronic devices 1 -. 
From a purely theoretical point of view, the dynamics of 
disordered magnetic systems represents a central prob- 
lem in non-equilibrium statistical mechanics. One of the 
central questions arising in the analysis of ferromagnetic 
systems, is the link between the domain structure and 
the hysteretic properties, such as the coercive field, the 
power losses and the noise. From the experimental point 
of view, it is possible to observe the domain structure 
on the surface of the sample by magneto-optical tech- 
niques, while the bulk behavior is accessible by induc- 
tive techniques. Several models have been developed 
to understand the experimental results, ranging from 
Ising type models considering the reversal of a set of in- 
teracting spins withSiMiSiS or without disorder^ ^ 10 ' 11 , 
to models focusing on the dynamics of a single domain 
walli^^iii^^ii 7 .^^. These two classes of models rep- 
resent two extreme situations: spin models are appro- 
priate when dipolar interactions are negligible. On the 
other hand, in soft magnetic materials magnctostatic ef- 
fects induce wide parallel domain walls that determine 
the magnetization properties. 

The dependence of the hysteresis loop area on the fre- 
quency and the amplitude of the applied magnetic field is 
often referred to as dynamic hysteresisi^ in the context 
of thin films and to power losses in bulk materials. In 
bulk materials, energy dissipation is dominated by eddy 
currents propagation 1 , while in thin films this effect is 
negligible. This led to the general belief that dynamic 
hysteresis in thin films is ruled by completely different 
laws than in bulk samples. The results of the models 
used to analyze dynamic hysteresis are often interpreted 



by assuming a universal scaling law for the dependence of 
hysteresis loop area A on the temperature of the system 
T, the applied field frequency to and amplitude Hq. The 
experimental estimates of the scaling exponents display a 
huge variability^^i^^^s^Ti^!., and the validity 

of a universal scaling law is still under debat e 21 ' 32 . In 
this respect, it was recently shown that the theory of loss 
separation developed for bulk materials could be equally 
well applied to thin films, since the precise nature of the 
damping (i.e. eddy currents or spin relaxation) does not 
change the basic equations 34 . 

The Barkhausen effect 36 consists in the irregularity of 
the magnetization variation while magnetizing a sample 
with a slowly varying external magnetic field. It is due 
to the jerky motion of the domain walls in a system with 
structural disorder and impurities 3 ^. Once the origin of 
Barkhausen noise (BN) was understood 3 -^, it was soon 
realized that it could be used as a non-destructive and 
non-invasive probe to investigate the magnetization dy- 
namics in magnetic materials. From a theoretical point 
of view, it is a good example of dynamical critical behav- 
ior, as evidenced by experimental observations of power 
law distributions for the statistics of the avalanche size 
and duration 3 -. There is a growing evidence that soft 
magnetic bulk materials can be grouped into different 
classes according to the scaling exponents values 3 ^. The 
Barkhausen noise is also an example of dynamics of a 
system presenting collective pinning when a quenched 
disorder is present, and it belongs to the family of the 
so-called crackling noises 3 . This kind of noise is exhib- 
ited by a wide variety of physical systems, from earth- 
quakes on faults to paper crumpling. So the relative ease 
to study crackling noise in magnets make the BN useful 
to get a deeper insight on different complex systems. 

The most successful models used to describe the 
BN have considered the dynamics of a single domain 
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walU^^^i&iL 1 ^, where the effect of the other do- 
main walls is accounted for by an average demagnetiz- 
ing held 16 . This simplification is justified because the 
Barkhausen noise is usually measured around the coer- 
cive held where the signal is stationary^. In this regime, 
the domain walls are reasonably distant from each other 
and one can assume that their mutual interaction is neg- 
ligible. The validity of this assumption, however, has 
never been established firmly: multi-wall effects may in 
principle affect the domain walls dynamics, and thus dy- 
namic hysteresis and BN. Even if we consider a single 
domain wall, its magnetostatic energy should depend on 
the number of the domain walls at least in an effective 
medium sense. In particular, the demagnetizing factor, 
that plays a crucial role in the BN statistics^, should be 
correctly evaluated only considering the entire domain 
structure. An attempt to consider the effect of wall-wall 
interactions has been made using a spring-block model 40 , 
but this model does not take correctly into account the 
long range nature of the dipolar interactions which un- 
derlie the multi-wall effects. 

In this article, we consider a system with many 
rigid parallel Bloch domain walls and study their mo- 
tion, driven by an external (triangular) magnetic held, 
in a disordered material with strong in-plane uniaxial 
anisotropy. The model is based on the interplay between 
the dipolar and the external held contributions, in the 
presence of structural disorder. Owing to the simplic- 
ity of the parallel Bloch conhguration, the magnetostatic 
energy can be treated analytically. We can calculate per- 
turbatively the demagnetizing factor function of 

the structural and geometrical parameters of the system. 
The perturbative calculation perfectly agrees with the 
results of the simulations thus representing a link be- 
tween a macroscopic measurable quantity and the mi- 
croscopic dynamics of the system. Moreover, dynamic 
hysteresis is investigated by integrating numerically the 
equations of motion of the domain walls. We analyze the 
behavior of the coercive held as a function of the applied 
held frequency, and of the disorder density and inten- 
sity. Furthermore, with our model we can investigate the 
Barkhausen noise. We hnd that the probability distribu- 
tions for the size, the duration and the amplitude of the 
Barkhausen avalanches show a power law behavior with 
a cutoff, in qualitative agreement with the available ex- 
perimental data. These results show that multi-domain 
effects are in principle important, since they give rise to a 
non-trivial BN and modify the effective parameters (e.g. 
the demagnetizing factor) describing the motion of a sin- 
gle domain wall. A more complete treatment should in- 
volve the dynamics of a system of flexible parallel domain 
walls, but this goal is beyond the scope of the present 
work. 

The manuscript is organized as follows. In section [TT1 
we present an overview on the energetics of films with 
parallel Bloch domain walls, and we compute magneto- 
static, disorder and external held contributions to the 
equations of motion of the domain walls. In Sec. IIIII 



we present the extended perturbative calculation of the 
magnetic susceptivity and the comparison with the sim- 
ulation results. Next, we present the results obtained by 
our simulations for the analysis of the dynamic hysteresis 
(Sec. [V]) and of the Barkhausen noise (Sec. IVI[) . Finally, 
our results are resumed in Sec. IVIII 



II. INTERACTIONS IN A MULTI-DOMAIN 
STRUCTURE 

Our purpose is to study the motion of n parallel Bloch 
domain walls in a disordered ferromagnetic system, un- 
der an external magnetic held driving. The aim is to 
write an equation of motion for each wall and integrate 
the system of equations numerically. To this end, wc 
calculate the total forces acting on the domain walls 
whose positions at time t are described by the vector 
x(i) = {xo(t),Xi(t),X2(t), ,x n {t)}, As we are inter- 
ested in the macroscopic response, we do not consider 
the details of the internal structure of the walls, and 
treat only the magnetostatic, the disorder and the ex- 
ternal held contributions. Thus the total force F(k,t) 
acting on the k-th wall at time t is given by: 

F(k,t) = F m (k,t) + F dia (k,t) + F ext (k,t). (1) 

In equation (JTJ) , the magnetostatic term F m (k, t) takes 
into account the interaction between the magnetization 
and the stray (magnetostatic) held generated by the mag- 
netic charges due to discontinuity of the magnetization 
vector at the upper and lower boundaries of the sample, 
Fdi s {k,t) models the contribution of structural disorder, 
impurities, defects and so on, and F ex t(k, t) describes the 
interaction between the magnetization and the external 
magnetic held. A detailed expression for Eq. |T]) is ob- 
tained by computing the energy E — E m + Edi S + E ext for 
a generic conhguration x and then deriving it according 
to: 



In the following subsections we will discuss these terms 
in more detail. 



A. Magnetostatic force 

We consider a sample with total length L, height 2d 
and thickness e, with n domain walls displaced in the po- 
sitions {xo(t),Xi(t), X2(t), , x n (t)} at time t, separat- 
ing Ti—l domains of alternating magnetization. We con- 
sider a system with strong uniaxial in-plane anisotropy 
along the z direction. The even domains have a magneti- 
zation equal to M s and the odd ones to —M s , being M s 
the saturation magnetization of the material (see Fig. [T] 
for a dehnition of the parameters). In order to calculate 
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the magnetostatic contribution F m (k, t) to the total force 
on the k-th domain wall at time t, we first compute the 
magnetostatic energy of the system E m for a generic ar- 
rangement of the walls and then derive it with respect to 




- (j. M s y2(-iy dx dy[cp(x,y,d) - <f>(x,y,-d)] 
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The scalar potential is given by the surface integral 



X „ X j X 2 



(j)(x,y, z) = J dS ^ 

n— 1 pxj + i re 

/ dx' / dy'M s 
j=o Jx j Jo 



y/(x - x') 2 + (y - y') 2 + (z - df 



\J {x - x 1 ) 2 + {y~ y') 2 + (z + df 



since the surface charge density a value is M s or — M s on 
the upper and lower interfaces (x, y, ±d) and elsewhere. 
The magnetization is uniform inside the domains, thus no 
volume density of charge is taken into account. 
Thus we obtain 



FIG. 1: (Color online) Sketch of the system parameters, for 
an array of wall positions {2:1,2:2, , x n }. 

The magnetostatic energy E m can be expressed in 
terms of the demagnetizing field as 

E m = -- fi j M H d dxdydz, (3) 

J sample 

where M is the magnetization vector and the integration 
is taken over the whole sample. Eq. (|3|) can be rewritten 

as 

n-l 



E m = -^o^(-l) 1 / dx I dy I dzM s (z-U d ), 



i=0 



Xi+l 



since the magnetization lies on the z axis. Using that 
Hd = — V^, where <j) is the scalar potential, it follows 
that 



E m = Mo M s 2 (- 1 

i,j=0 



fX i + 1 re re 

1 1 I dx dx' dy dy' 

Jxj JO JO 



yj(x - x') 2 + (y- y') 2 yj{x - x') 2 + (y - y') 2 + 4d 2 

(4) 

Assuming that e << 2d, as in most of the samples of 
experimental interest, we can neglect the term (y — y') 2 
with respect to 4d 2 in the Eq. [4j Therefore, the result of 
the integral is given by 



n — 1 ^ 

E m = ^M* £(-ir + ^ 



i,3=0 



z-Hd = -T^(f>(x,y,z), 



and therefore 



{[g(x i+1 ,x j+1 ) +g(xi, xj) - g(x i+1 ,Xj) - g(xi,x j+1 )] 



n-l 



i=0 



E m = 5Mo2_/( _1 ) / dx I d y I dzM s —ct>(x,y,z) 







[f(x i+ i,X j+ i) + f(Xi,Xj) - f(x i+ i,Xj) - f(Xi,X j+ i)]} 



(5) 



where 
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[ f(x,y) = f(y,x) = 2^/ Ad 2 + (x - y) 2 



y) in 



y — x+^J Ad 2 + (x — y) ; 



x-y+\f id 2 + (x-y)' 



g(x, y) = g(y, x) = {x- y)e 2 In 



y-x+y/e 2 + (x-y) 2 



-§(o> 



y i3 + (x _ )2 eln / V(fz^!±!!- e N ) 
-y) 2 - 



,3/2 



To give an idea of the behavior of the magnetostatic 
energy, in Fig. ((2|) we have plotted E m (t) for a sample 
of fixed length L and d and for three values of the thick- 
ness, e = 0.1; 0.01; 0.001, as a function of the number 
of (periodic) domains. As it can be seen, from a purely 
magnetostatic point of view, the energy is a decreasing 
function of the number of domains. In fact, the total 
number of domains in a sample is determined by the in- 
terplay of the magnetostatic term and the domain wall 
energy, linearly increasing as a function of n — . 
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FIG. 2: (Color online) Magnetostatic energy E m as a function 
of the number of domains n, for three different values of the 
system thickness e = 0.1; 0.01; 0.001, and for L — 1., and 
d = 10, calculated for even numbers of domains. Solid lines 
are guides to the eye. 

The magnetostatic force on the fc-th wall at time t is 
thus 



F m (k,t) 



dE m 
dx k (t) 



n-l 

2 Al0 A/ s 2 £(-l) 4+fc 



dg(x i+1 ,x k ) dg(xi,x k ) 



dx k 



dx k 



df(x i+ i,Xk) df(xi,Xk) 



dx k 



dx k 



(6) 



df(x,y) _ df(x,y) = [ y-x + ^/4d 2 + (x - iff 



dy 



dx 



-y+ yj±d 2 + {x - yf 



^§^=-^^=2( X -yWe 2 + (x-yy 



+2e(x — y) In 



y/e 2 + (x -y) 2 -, 

V e2 + ( x -y) 2 + ' 



; ln 



y — x + \A 2 + ( x - yf 



x-y+ y/e 2 + (x - y) 2 



-2(x-y) 2 sign(x~y). 

(7) 




where 



FIG. 3: (Color online) Magnetostatic force F m (k,t) on the 
fc-th domain wall as a function of fc (red dots), for a periodic 
configuration of the domain walls positions. The sample has 
n = 21 domains, e = 0.002, d — 10, and L = 1. The solid line 
is a guide to the eye. 

In Fig. we show the magnetostatic force F m (k,t) 
calculated at each wall k for a periodic configuration of 
n = 21 domains. As it could be seen, walls with k of 
opposite parities are driven in opposite directions, and 
the absolute value of F m (k,t) is higher for the external 
walls, which is a finite size effect tending to reduce the 
size of the boundary domains. For a large odd number 
of domains, in the center of the sample, all the F m (k,t) 
would be equal in magnitude, as it can be seen in Fig. 
Q. For an even number of domains, the central wall 
experiences a zero force for symmetry reasons. As ex- 
pected, for a generic (even or odd) number of domains, 
the effect of the magnetostatic force is to move the walls 
in the direction that decreases the magnetization, in or- 
der to minimize the magnetostatic energy [see Eq. ([5])]. 



B. Disorder 

Different sources of inhomogeneities are found in virtu- 
ally all ferromagnetic materials. The presence of struc- 
tural disorder is essential to understand the hysteretic 
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behavior, and especially to account for the residual co- 
ercive field when the frequency of the external driving 
field vanishes. Disorder is provided by vacancies, non- 
magnetic impurities, dislocations or grain boundaries in 
crystalline systems, variations of the easy axis between 
different grains for polycrystalline materials, and internal 
stresses for amorphous alloys. 

We consider here only quenched (frozen) disorder that 
does not evolve on the timescale of the magnetization 
reversal, which is usually a realistic approximation for 
ferromagnetic systems. At the beginning of a simulation, 
we extract the positions of a fixed number n p of pinning 
centers with a uniform probability distribution all over 
the sample. Each pinning center interacts with every 
domain wall by a located potential of the type 

U p (x) = Aexp[-(x/0 2 ], (8) 

where x is the distance between the pinning center and 
the wall, £ is a correlation length and the amplitude A is 
extracted from a uniform distribution, considering that 
the strength of the pinning centers should vary for struc- 
tural reasons. In this way we produce a disordered land- 
scape in which the domain walls evolve. 

C. External field 

The interaction between a ferromagnetic system and 
an external magnetic field H is described by the energy 
term: 

E ext = -^ [ M • H dV. (9) 

J sample 

In our model the uniaxial anisotropy and the external 
time dependent magnetic field H ext (t) are parallel to the 
easy axis of the sample, and thus parallel or anti-parallel 
to the magnetization. Thus in this case the total energy 
is 

71-1 

E ext = 2defi M s H ext (x t+1 (t) - Xi(t)) , 

i=0 

and the external field contribution F ext (k, t) to the force 
F(k,t) on the fc-th wall at time t is given by the partial 
derivative of Eq. flSJ) with respect to the wall position 
Xfc(i), and could be thus written as 

F ext (k,t) = 4ti M s deH ext (t)(-l) k+ \ (10) 
which depends of course on the parity of the domain wall. 

III. DEMAGNETIZING FACTOR AND 
MAGNETOSTATIC SUSCEPTIVITY 

An interesting analytical result that could be obtained 
from our model is the dependence of the demagnetizing 



factor k on the number n of the parallel domain walls. 
The parameter k enters as a mean field term in several do- 
main wall models used to describe the B N 13 ' 14 i 15 ' 16 i 18 . Its 
precise dependence on the sample geometry is usually ap- 
proximated considering a uniformly magnetized sample 
and the domain structure is ignored. Understanding this 
point is important since k determines for instance the size 
dependence of the Barkhausen avalanche distribution 39 . 

The demagnetizing factor k can be rigorously defined 
for uniformly magnetized ellipsoid from the relation = 
— kM. In more general situations, we can define it as an 
effective quantity linking the average demagnetizing field 
to the magnetization k = —{Hd)/M. Here we compute 
k from the magnetic response of the system to a small 
perturbative external field. Since at equilibrium a small 
change in the external field is compensated by an increase 
in the demagnetizing field (Hd + H ext = 0), we can link k 
to the magnetostatic susceptivity x = dM/dH ext by the 
relation k = 1/x- 

For an extended system with a large and even number 
of domains n, the equilibrium configuration consists of a 
domain wall array 

Xi = — Vt 
n 

(note: in what follows we will omit the t dependence of 
the functions). We can choose a positive perturbative 
field SH ext and thus consider the perturbed walls posi- 
tions: 





iL 


Xi 






n 


x 


= 


•i' 11 


= L 



that correspond to a (small) increase of the magnetiza- 
tion dM ~ nu, with u > 0, uniformly distributed over 
the domain walls. Since the calculation of the magnetic 
susceptivity is quite involved, we discuss it in the Ap- 
pendix and report here only the final result: 



dH ext k eA(n) ' 

where 

A(n) = 2 [7 + ln(n)] , 

and 7 ~ 0.577215. Higher order corrections in n have 
been neglected. 

In Fig. 2] we compare the results of our simulations 
and the theoretical result of Eq. [TI] for a system with 
d = 400, e = 0.01 and L — 1, without disorder, by study- 
ing k as a function of the number of domain walls n (coin- 
ciding with the number of domains in the limit n — > 00). 
The simulations are performed starting from a periodic 
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configuration of a high and even number of domains, and 
letting the walls relax under an external field that is very 
low, in order to have a displacement from the equilib- 
rium position of the walls of the order of some per cents 
of the distance between nearest-neighbor walls. To fulfill 
the hypothesis under which Eq. [ITJ is calculated, we use 
here periodic boundary conditions. As it can be seen, 
the agreement between theory and simulations is excel- 
lent, even if the simulations have been performed for fi- 
nite samples. The value of the magnetic susceptivity is 
of the order 1/k = x~10 4 , which is a realistic value for 
a ferromagnetic system. It would be very interesting to 
test the prediction of Eq. [Tl] with experimental results 
on materials with low disorder. Finally, it is interesting 
to notice that the demagnetizing factor k = 1/x displays 
an inverse dependence on the sample length d and a lin- 
ear dependence on the thickness e. Hence, as expected, 
reducing the sample aspect ratio leads to an increase in 
k. This fact has been exploited in Ref. [3^ to tune K and 
study its effect on the BN. 

2.5x10"' 



2.0x10" 
K 



1.5x10 



o simulations 
— K=EA(n)/d 



Kutta algorithm. The external magnetic field is a saw- 
tooth signal with rate u>, being H cx uj. 

We always start our simulations with an M = at 
H ex t — periodical configuration, and drive the sample 
to positive and then to negative saturation. We include 
the possibility of nucleation and annihilation of the do- 
main walls by fixing a minimum interaction range be- 
tween two nearest-neighbor walls, namely <5 TO , n . If two 
walls get closer than <5 mi „ , we stop them and consider the 
pair of walls as annihilated. This means that their con- 
tribution on the sum involved in the calculation of long- 
range magnetostatic force (see Eq. [5]), and thus their 
effect on the equations of motion of the other walls, is no 
more considered. Besides, we go on calculating the total 
force on the annihilated pair of walls. When this force be- 
come strong enough (and of the correct sign) to let them 
escape from the range 5 m i n , we nucleate them and con- 
sider again their contribution to the magnetostatic force 
calculation, Eq. [S] This corresponds to have a nucle- 
ation barrier equal to zero. Since we are interested in the 
study of Barkhausen noise and of the coercive field be- 
havior connected to dynamic hysteresis, this choice does 
not influence our results, as we have checked. In fact, 
all the phenomena that we are interested in take place 
in the central region of the hysteresis loop, namely the 
one close to the coercive field H c , while the choice of the 
nucleation barrier involves only the region of the cycle 
close to saturation. Moreover, for the same reason, our 
results are independent on the choice of S m i n , as long as 
it is chosen in a reasonable range. 



COERCIVE FIELD AND DYNAMIC 
HYSTERESIS 
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A. Dynamic hysteresis 



FIG. 4: Demagnetizing factor Ki ELS cL function of the number of 
domain walls n (semilogarithmic scale): comparison between 
simulation results and perturbative calculation (dashed line) 
for d = 400, e — 0.01, and for unitary L. 



IV. NUMERICAL MODEL 

In the previous sections we have calculated the various 
contributions to the total force acting on each wall (Eq. 
[JJ thus it is now possible to simulate the dynamics of 
the system. Since the number of the walls is n, and 
we consider in this part closed boundary conditions (the 
boundary walls do not move), we have to write down n— 2 
equations of motion given by 

dx i* 

T-^-=F(k,t)=F m (k,t)+F ext (k,t)+F dis (k,t). (12) 

In the following we set the damping coefficient T to unity. 
We integrate these equations using a fourth order Runge- 



We first consider the effect of the external field rate 
on the hysteretic behavior. In Fig. E{a) we show the 
hysteresis loops obtained for various rates for a system 
with n = 20 domains (since from now we will deal with 
a reasonably high number of domains, we will use with 
a slight abuse of terms n for the number of domains), 
where m is the normalized magnetization m = M/M s 
and h is the normalized magnetic field h = H/M s , As 
expected from experiments and general considerations, 
small (high) frequencies correspond to narrow (large) cy- 
cles. To quantify this observation we can focus on the 
coercive field H c behavior. In Fig. [5fb) we show the 
dependence of H c on the field rate u>. 

For all the considered values of n in Fig. [SJ n = 
20,40,60 (for a discussion of the behavior of H c as a 
function of the number of domains n see the next sub- 
section), H c shows an increasing dependence on u>. The 
H c dependence on lo is quite close to a power law of the 
form 

H c = H p + Auj 1/2 , (13) 
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0.5 - 



E 



-0.2 



0,035 




n 


20 


40 


60 


H p 


0.024 


0.019 


0.016 


Cex 


0.003 


0.003 


0.0036 


r 


6622 


4752 


2989 


A 4 


183 


153 


173 


n 


18 


21 


32 


Vb 


0.0003 


0.0003 


0.0002 



TABLE I: Fit parameter values for Fig. [5jb) . 



interpret as the pinning dominated quasi-static compo- 
nent due to structural disorder, while the behavior of 
H c in the high rate regime represents the domain wall 
dominated dynamic contribution. Moreover, the scaling 
function of Eq. [14] describes the evolution from the adia- 
batic to the domain wall dominated power-law behavior 
without invoking mechanism crossover between domain 
wall propagation and nucleation of new domains. In fact, 
even if the nucleation process is included in our model, 
besides in a very simplified form, nucleation takes place 
only in the regions of the hysteresis loops very close to 
saturation. Thus the cycles are very oblique [see Fig. 
Ufa)] . Such a kind of hysteresis loops are expected to be 
dominated by domain wall motion. 



5x10 

dh/dt 



FIG. 5: (Color online) (a) Hysteresis loops for different ex- 
ternal field rates, for a system with n — 20 domains, and (b) 
Coercive field as a function of the external field rate for three 
different system sizes, n = 20,40,60; h c vs dh/dt in a linear 
plot, in comparison with a square root fit (dashed line) and 
the fit with the power loss formula (continuous line, see text). 
Every point is mediated over 100 realizations corresponding 
to different disorder configurations. 

as suggested for experiments on permalloy thin films and 
microstructures fabricated from the same sample^, and 
by the theory presented in Ref. 112 . 

Nevertheless, as it can be seen in FigJSfb), our curves 
are best fitted by the loss separation formulaic, 

H c = H p + C ex [{l + rH) 1 / 2 - l], (14) 

where C ex = noVb/2, and r = 4r/z/(ngVb). Here no is the 
number of active walls in the quasi-static limit (to — > 0), 
Vb is a characteristic field which controls the increase of 
no due to the excess field, fi is the permeability and H p 
is the static (hysteretic) component. The fit parameters 
values are shown in table HI 

The behavior of H c as a function of u), suggested in 
Eq. 1141 means that in the adiabatic limit (low frequen- 
cies) H c goes to a non-vanishing value H p that we can 



B. Role of the number of domains 

Another interesting issue to analyze is the effect of the 
number of domains n on the system behavior. We re- 
mind here that at the beginning of each simulation we 
set the maximum number of domains of the specimen. 
This number can decrease down to one while driving the 
system to saturation, and then, increasing the external 
field, it is completely restored around the coercive field. 
From a qualitative point of view, keeping constant the 
other system parameters (i.e. the physical dimensions of 
the specimen under study and the pinning centers density 
and strength /□), a high (low) number of domains cor- 
responds to hysteresis loops with a lower (higher) per- 
meability, in accord with Eq llll as it can be seen in 
Fig. E)Ja). This happens because the increase of the 
domain number increases the relative relevance, in the 
total force balance, of the dipolar contribution. In fact, 
an absence of collective behavior leads to square shaped 
cycles, where the whole magnetization reversal process 
takes place in a few avalanches. Otherwise, in materials 
in which dipolar interactions are relevant, oblique hys- 
teresis loops are expected. In our simulations, the coer- 
cive field H c decreases as a function of n, as expected [see 
Fig. Mb)}. 
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FIG. 6: (Color online) (a) Some hysteresis loops for different 
system sizes, n — 20, 40, 60, for an applied field rate u — 
1/20, and (b) Coercive field H c as a function of the number 
of domains n. Every point is an average on 100 samples, 
corresponding to different disorder realization, for an applied 
field rate uj = 1/1000. The line is a guide to the eye. 



C. Effect of disorder 

As we have already mentioned, the presence of struc- 
tural disorder has a crucial role in the physics of ferro- 
magnetic systems. Disorder pins the domain walls, that 
thus are not able to move until the external field reaches 
a sufficient value to overcome the pinning and let the 
wall move again. Hence we will expect, from this simple 
argument, that more disordered systems have larger hys- 
teresis loops compared to purer specimens, for the same 

UJ. 

We remind here that it is difficult to quantify the dis- 
order contribution in real systems. In our model we have 
two degrees of freedom to control the disorder: we can 
tune the density of pinning centers (by controlling the 
total number of pinning sites n p ) and their strength /o. 

In Fig. [7] we show some hysteresis cycles for varying 



FIG. 7: (Color online) hysteresis loops for varying: (a) disor- 
der strength /o (n p = 600) and (b) number of pinning centers 
n p (f = 1.2). The rate is u = 1/50. 



values of the two parameters /o [see panel (a)] and n p 
[see panel (b)]. As we can see, their effect on the loops 
is qualitatively similar, and the behavior for increasing 
disorder (i.e. larger cycles for larger disorder) is in agree- 
ment with the general argument explained above. 

In Fig. [8l we show the coercive field behavior for in- 
creasing disorder, intended here as increasing both the 
values of n p and /o : as expected, stronger disorder im- 
plies higher coercive field. Moreover, it is important to 
notice that even varying the disorder strength or density, 
the H c dependence on the square root of the applied field 
rate does not change. 



VI. BARKHAUSEN NOISE 

The Barkhausen noise is an unavoidable feature of 
magnetic hysteresis in disordered samples and its statis- 
tical properties have been widely studied in experiments 
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FIG. 8: (Color online) The reduced coercive field h c depen- 
dence on the external field rate dh/dt for a system with n — 40 
domains and two different number of pinning centers n p and 
disorder amplitude /o. The applied field rate is ui — 1/1000. 
Data are fitted by Eq |14l Fit parameters are 
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FIG. 9: (Color online) Time sequences of Barkhausen noise. 
The labels indicate the applied field rate. We report portions 
of simulated signals corresponding to the same time interval. 
The time is measured considering an integrating time step 
At = 0.05. 



and models 3 - 9 .. The experiments in soft ferromagnetic 
materials have been divided into universality classes de- 
pending on the avalanche distributions 39 . Experimental 
results are well described by models involving a single 
flexible domain wall, and ignoring any interactions be- 
tween different walls. Different universality classes are 
set by the dominant interactions between parts of the 
walls. Here we explicitely ignore these local interactions 
and want to account only for the interactions between 
walls. 

In Fig. [HI we show some simulated signals for various 
rates of the external field. Partial regions are chosen for 
the signals, corresponding to the same time interval. As 
it can be seen, our results reproduce a well-known charac- 
teristic of the experimental Barkhausen noise: at low fre- 
quencies, the signal is composed of separated avalanches, 
while this property is progressively lost at higher frequen- 
cies. In the latter case, the signal appears as a continuous 
sequence of peaks and the single avalanches are no more 
distinguishable. We have analyzed the Barkhausen sig- 
nal considering the distribution of the signal amplitude, 
avalanche sizes and durations. 

The amplitude probability distribution of the magnetic 
induction flux rate <f> is a measurable quantity and thus 
has been often used in order to test the reliability of a 
model for Barkhausen noise (see for instance Ref. Il31). In 
fact, the flux rate is related to the velocity of the domain 
walls v. In Fig. fTUf a) we show the rate dependence of the 
probability amplitude distribution. In the non-adiabatic 
limit, in fact, the domain walls velocity is expected to 
depend on u>. We notice that P(v) passes from a almost 
symmetric shape for very high frequencies to an asym- 



metric one for low frequencies. 

In the ABBM modeli 3 -, the predicted shape of the am- 
plitude probability distribution is a Gamma function 

p^airv™^. (15) 

We remind here that this formula is obtained under the 
assumption that the disorder is a correlated (Brownian) 
process. 

The mechanism giving rise to the Barkhausen effect is 
the domain walls motion, and that a domain wall motion 
can be described as a stationary Markov process, associ- 
ated with fixed and well-defined values of the magnetiza- 
tion rate imposed to the system and of the permeability 
of the specimen fi, associated with the part of the hys- 
teresis loop where the domain wall motion is assumed to 
take place. Even if our model is different from the single 
wall model ABBM, the formula of Eq. Qj)]is able to fit our 
probability distributions P(v) for high rates, although it 
does not supply a suitable explanation for the low rate 
regime [see Fig. fTOT b)]. 

Next, we evaluate the probability distribution for the 
duration T and the size S = J T dtv{t) of the avalanches. 
The distributions are shown in Figs. [TTTa) and fTTTb). 

As it can be seen, both distributions display a power 
law behavior of the type P(S) ~ S~ T and P(T) ~ T~ a , 
for almost two decades for the maximum studied size 
(n = 60), with a cutoff. The scaling exponents of these 
distributions are respectively t ~ 1.1 and a ~ 1.2. These 
exponent values, however, have never been reported for 
bulk soft magnetic materials where parallel domain walls 
are commonly observed 3 - 9 -. This is not surprising as we 
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FIG. 10: (Color online) (a) Probability distribution of the 
signal amplitude (i.e. of the domain walls velocity) v, for 
various frequencies u of the external magnetic field (u = 
1/1300,1/700,1/400,1/300,1/200,1/100, for a system with 
n = 60 domains, (b) Probability distribution of the signal 
amplitude for n = 60 and lo = 1/130, fitted with the ABBM 
function, P(v) — av c ~ 1 exp(—cv/ <v>). 



neglect local interactions, and any deformation of the 
domain wall. It is nevertheless remarkable that just 
the combined effect of dipolar interactions and disor- 
der yields power law avalanche distribution in a multi- 
domain model. The latter appear, however, essential to 
recover quantitatively the experimental results in bulk 
samples 39 . To completely resolve this issue, it would be 
necessary to study the dynamics of an array of flexible 
and interacting domain walls. 



VII. CONCLUSIONS 

The technological and theoretical relevance of the hys- 
teresis in ferromagnetic systems motivates the effort in 
modeling the phenomenon. Of particular interest is the 
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FIG. 11: (Color online) (a) Probability distribution of the 
avalanche size S, and (b) Probability distribution of the 
avalanche duration T, for three different system sizes (n = 
20,40,60). 



study of the dynamic hysteresis, connected with the 
power losses of the material, and of the Barkhausen noise, 
which represents a tool for the investigation of the mag- 
netization dynamics. 

In this article we have presented a model for the 
dynamics of a system of parallel and rigid uncharged 
(Bloch) domain walls, moving in a disordered landscape 
driven by an external magnetic field. The analyzed con- 
figuration of domain walls is a common case both in two 
and three dimensional samples. Our model is based on 
the interplay between the long-ranged dipolar interac- 
tions and the external field and disorder contributions. 
Due to the simplicity of the wall configuration, we cal- 
culated perturbatively the demagnetizing factor and the 
magnetostatic susceptibility as functions of the geomet- 
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rical parameters of the system, namely the thickness, the 
height and the number of domains at equilibrium, in the 
absence of external magnetic field. The simulation results 
perfectly agree with the calculations, and we have thus 
been able to link a measurable macroscopic quantity to 
the microscopical system behavior. Next, we have inves- 
tigated the dynamic hysteresis by means of the study of 
the coercive field behavior, obtained by simulating the 
equations of motion of the walls array. The coercive 
field displays a power law behavior as a function of the 
applied field rate, with an exponent in agreement with 
experimental available data. Finally we have analyzed 



the Barkhausen noise, constructing the probability dis- 
tributions of the Barkhausen avalanches size and time 
duration, and found power laws with cutoffs for both the 
distributions. Moreover we find that the probability dis- 
tribution of domain walls velocity deviates at low driving 
rates from the ABBM related formula. Since the latter 
model is based on the motion of a single domain wall, we 
can conclude that the multi-domain effects are important 
in ferromagnetic systems, since they affect system param- 
eters like the demagnetizing factor and the behavior of 
the Barkhausen noise. 
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VIII. APPENDIX 

We can develop the magnetostatic force F m (k) on a 
generic wall k (see Eq. [6]), around its equilibrium posi- 
tion: 

F m (k) = F^{k) + 5F m {k) = 6F m (k), 
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since the equilibrium term F^(k) is by definition. It 
follows: 



Fm(k) = 

3=0 



dF m (k) 



Uj—0 



n-l 

E 



dF m (k) 



(-iy +1 uj, 



(16) 



where the last equality descends from the fact that the 
boundaries of the sample are fixed. 

For simplicity in the calculations, let us separate the 
contribution k = j from the others, and rewrite Eq. [TB] 
as: 



f d 2 f(x,y) 
dx 2 



V4d 2 + (x - y y 



d 2 g{x,y) 



dx 2 



-2e In 



(Ve 2 + (x- y) 2 ~\x- yfj (18) 



^e 2 + (x - y) 2 - e 
y / e 2 + (x-y) 2 + e 



In order to handle analytically Eq. [TTJ let us restrict 
our study to the limit e << \x — y\ << L << d, which 
is the physical case in most ferromagnetic samples. The 
second derivatives of Eq. [18] thus become 



F m {k) = — u k { 

ox k 



-i) fc+1 +E 



dF m (k) 
dx 3 



M,=0 



2M^ (-l) k+1 u k 



n-l 



i=l 



d 2 g(xj+i,x k ) _ ^ 2 d 2 g(x j ,x k ] 



dx 2 



dx 2 



2 i'9 2 f{x j+1 ,x k ) 2 d 2 f(xj,x k ) 



dx 2 



' d 2 f{x,y) 
dx 2 

d 2 g(x,y) 



1 

"d 



2r 2 



9a; 2 | a; — y\ 

and F m (k) could be in this limit rewritten: 

2e 2 



F m (k)=Afi M 2 J2 



(19) 



+ 4Af 2 g(-D J+fc f^%^- e 2 



3=0 



dx 2 



x(-iy +k ((-i) k u k -(-iy Uj ) 



(20) 



where we have used the symmetry properties of the 
derivatives (see Eq. [7]). In the last sum we have included 
again the fc-th term since it is null. If we neglect the con- 
tribution from the boundaries, that could be done since 
we are in the n — ► oo case, we can translate j ; + 1 — > j in 
the first sum, together with the sum index. We obtain: 



To have some insight on Eq. [20l we can study the effect 
of the magnetostatic force in two simplified situations: 

• Let us imagine that all the walls are fixed but the 
fc-th. Thus Eq. [20] can be rewritten: 



F m (k) = 4ii M 2 u k J2(-iy 



d 2 g(xj,x k ) 2 d 2 f(x J ,x k ) 



dx 2 



dx 2 



F m (k) = -8 M oM s 2 e 2 



(-1)' 



Xk\ 



r(-i) fc «fc 



\%i %k | 



-A^M s {-iy^{~iy U] 

3 

So finally, it results 



d 2 g(xj,x k ) 2 d 2 f(x J ,x k ) 



dx 2 



dx 2 



F m (k) = 4/i M. 



n-l 

' 2 E 

3=0 



d 2 g{xj,x k ) 2 d 2 f(xj,x k ) 



dx\ 



dx 2 



x ((-iyu k -(-i) k Uj ), (17) 

where the double derivatives of Eq. [7] are: 



that means that every wall i contributes with a 
term of sign (— l) l+1 . Since the most important 
contribution in the sum of Eq. [201 is due to the 
nearest neighbors, if k is odd (even) and has thus 
a positive (negative) perturbation on the position, 
it will be subject to a negative (positive) magne- 
tostatic force, recalling the wall to the equilibrium 
position. 

• Considering the case with all the perturbations 
Ui = u, Vz. The magnetostatic force becomes 



F m (k) = 4M 2 e 2 



x k \ 



d, 
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x (-iy+ k ((-i) fe - 



(21) 



So only the walls with the opposite parity with 
respect to k will contribute to F m (k). Thus the 
even (odd) walls, that have a negative (positive) 
displacement, will experience a positive (negative) 
force, recalling the wall to the equilibrium position. 

Let us restrict ourselves in the following to the limit 
above, Ui = u, Vz. So, starting from Eq. [21] we can write 



F m (k) 



4fi M^€ 2 u ^ 



((-ir-(-i) fc ) 



by using the sum of an harmonic finite series of odd 
terms, where 7 = 0.577215. So finally 



F 



l (k)=^ 0A f 2un (~l) k A(n) 



where 



A(n) = 2 [7 + In + 2 In (2)] = 2 [7 + ln(n)] . 
The energy term due to the external magnetic field is 
given by Eq. [9l and in this case could be written 

Eext = -4n M s H ext de ^ Uj. 



i 

using the definition x\ — iL/n. Since L << d and n/\i ■ 
k\ > 1 V i, k, we can approximate 

Wl-^E^K-'M- 1 )') 



Therefore, F ext (k) is given by: 

dE ext 4fi M s H ext de £) 4 Ui 



F ext {k) 



dxh 



(-l) k+1 u k 



= (-l) K+L 4fi M s H ext de. 

Since at the equilibrium the total force F = 0, in absence 
of disorder it must be 

F m (k) + F ext (k) = 



defining a — S^lqM 2 e 2 un/L. In the sum 

*m(*) = a(-l)*E|^ ((-!)*"-!) 



the only non null terms are the ith terms i' with opposite 
parity with respect to k, that lead to a contribution 



F m (fc) = a (-l) fe (-2)£^ 



k\ 



Since n — > 00, we could consider k in the middle of the 
sample and write 

n/2 n/4 

Fm (k) = ~M~i) k E - = -M-i) fc E — x 
¥=1 t=i 



2a(-l) k [7 + In Q +2 In (2) 



(-l)^(n) + (-l) fc+1 4 Mo M s ff ext de, 



from which follows 



2euA(n)M s n 
Id ' 



The magnetization M could be written 
AudtM^n 2uM,n 

M = 



2dLe L 
So finally we can calculate the magnetic susceptivity 



dM M d 



X 



dH ext H ext eA{n) 
which was our goal. 



(22) 



